INITIATION OF MOTION : 

EFFECT OF ROUGHNESS CONCENTRATION 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 
MASTER OF TECHNOLOGY 


By 

MOHD. ASLAM 


to the 

DEPART MlENT OF CIVIL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

may, 1381 



II.T. 


k I f - 'j, « 

vL V. 


life, i 

•6 MAYfSjj 


;. 65977 


CB^lW 


IHI 



i 


GERTIFIGATE 


The present work entitled 'INITIATION OF MOTION; 
EFFECT OF ROUGHNESS OONCENTRATION ' by Mohd, Aslam is 
hereby approved as a creditable report on research carried 
out and presented in a manner which warrants its acceptance 
as a prerequisite for the degree of MASTER OP TBCHNOIOGI. 
The work has been carried out under my supervision and 
has not been submitted elsewhere for award of a degree. 


May 2, 1981 


d* 


J\0 V C‘-l a 


/u 


( Dr. T, G-angadharaiah ) 
Assistant Professor 
Department of Civil Engineering 
Indian Institute of Technology , Kanpur 



ii 


kamom-spmMim 

I wisli "to express my sinceres'fc and deep sense of 
gratitude to my thesis supervisor Dr, D. Gangadharaiaii , 
for his valuable guidance, constant inspiration and 
encouragement during the course of the present work. 

Sincere thanks are extended to Dr,. V. Lakshminarayana, 
Dr. K, Subramanya and Dr. S. Surya Rao of the Department 
of Civil Engineering for the thought provoking courses me 
received from them during the degree programs. 

I acknowledge with thanks to Mr. D.D, Porey for 
his inspiration and keen interest, which made the completion 
of the present study. 

Technical help provided by staff of Hydraulics 
and V/ater Resources Laboratory of the Department of Civil 
Engineering is greatfully acknowledged. 

And finally, the unreserved help received from 
Sri S.P, Sharma is also acknowledged. 


Mohd. Aslam 



ABSTRACT 


Tile flow conditions corresponding to the initiation of 
motion Of sediment grains in a bed having densely packed with 
ooheeionless sediments have been investigated very systamatically 
in the paet , Threshold flow conditions corresponding to a 
particular grain size, resting on a particular level with 
respect to general surface of the bed in a less dense bed 
have not been investigated. An attempt has been made to 
investigate for these flow conditions in the present work. 

The experimental programme of the present work 
covers three aspects j namely, the effect of roughness concen- 
tration, effect of size variation of a particular grain 
incomparision to grains constituting bed roughness and the 
effect of relative position of the grain of same size that 
of roughness elements on their initiation of motion. The 
following are the significant outcomes of this investigation. 

The drag co-efficient calculated from the mean v elo city 
data found to increase with decrease in roughness concen- 
tration, X t and tends to be independent of roughness 
concertration for X<0.1 , The theoretical bed level found 
to coincide with the average surface level of the rough bed. 

The decrease in roughness concentration, the increase 
in grain size, and the increase in elevation of grain, 
result in the reduction in critical shear stress, \ based 
on the average bed shear stress. However, critical shear 
stress e xpressed interms of shear stress, offered by 

a particular grain under investigation agrees with Shields 
critical shear stress value. 
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Constants 

Constant 

Diameter of glass bead 

Center to center spacing of glass bead. 

Variable diameter of glass bead. 

Acceleration due to gravity 

Average h. f'htcf roughness element 

Nikuradsc’* equivalent sand grain roughness 

V^d Shear Kanolds number, 

V 

Suffix s for rough, bed 
Suffices for smooth bed 
Average velocity atany height 
Shear velocity = v'Ty'p 

Intercept of law of wall 

Ordinate, distance measured perpendicular to 
the bed from the top of grain. 


Y^= i.Y'f ej Perpendiculer height from theoritical bed level. 



Greek alphabets 


5 Boundary layer thickness measured from the top 

of the grain. 

e The location of the apparent origin from velocity 

distribution above the smooth surface of the flat 
bed. 

Von Karman*s universal constant 
p Mass density of the fluid 

Pg Mass density of glass bead 

shear stress 

Hon- dimen sionalised critical shear stress. 

xjf Non-dimensionalised critical shear stress 

offered by single particle 

Kineiaatic viscosity 

^ Roughness concentration 



CHAPTER - 1 


INTRODUCTION MB LITERATURE REVIEW 
'••'I Introduction : 

Initiation of motion of sediment is a starting 
step of sediment transportation. An extensive work Has 
been carried out by various investigators to find the 
general conditions at #iich. the whole sediment bed is 
subjected to motion. In cohesionless bed, the sediment 
motion occurs by the movement of individual sediment 
grains. The knowledge of flow characteristics on the 
initiation of motion of grains in a densely packed sediment 
bed are available due to a systematic investigations 
carried out in the past. It is not known how the flow 
conditions at the initiation of motion of sediment 
particle in a less dense bed differs from that densly 
packed, bed, Further, very little information is 
available on how the relative size and position of grain 
alters the critical condition, A knowledge on these 
characteristics are in need for proper understanding 
of the mechanics sediment transportation. 
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The fluid flowing over a bed of cohesionlfiss 
sediment, exerts a fluid dynamic force on the grains and 
tends to move or entrain them. The force resisting the 
motion of Ci.^hesionless sedimoit is mainly due to their 
weight. When the fluid dynamic forces ^Drag force and 
uift forced acting on a sediment grain attains a particular 
vaLie and if that value is increased slightly, the grain will 
be subjected to the motion. This state of motion of sediment 
is called 'Initiation of motion’ ' Criti«al conditin’ or 
’threshold condition' . This threshold condition is a 
result of combined action of drag h lift forces on the 
submerged of the grain in cohesionless sediment bed. 

Tile dr;-, and lift forces acting on the grain depends upon , 
the roughness charactiristics of the bed and size of grain 
and its relative position with respect to. the average bed 
surface.. The relative spacing of the grains on the bed may- 
be represented as roughness concentration, which is 
defined as the ratio, of the projected area ofv raugtoess 
elements perpendicula.r.- to the flow to the floor area* 



In /cstigations on ^nltist-lun. of uotribn wre gene;-- 
nally carriedout fot the condition ,. -5 of the general motion 
of grains in the densely packed sediment bed. A few 
attempts have been made to investigate the threshold 
condition of a single grain resting either on densely packed 
sediment grain bed or on the smooth bedi Almost no invest- 
gation has been carried out to find the influence of 
rouglTn..'SS concentration and size variation on the threshold 
conditii.-n of an individual grains. In the present invest- 
igation an attempt has been m de to study the influence 
of roughness concentration, size v;nriation and its relative 
position on the threshold condition of the grain under 
observation. 

, 2 .^itetature Jiev i . 

Initiation of motion is one of the major aspects of 
sediment transport. This problem of initiation of motion 
has been analysed by several investigators and attempts 
have oeen mcode to develop relation between critical tractive 
force, sediment size and sediraent characteristics such as 
specific weight. The best known and widely accepted theory 
on the initiation of motion was proposed by Shields \1936;. 
Vanoni v1977; summarizes the recent contributions of various 
investigators on this aspect of the problem. A brief review 
of litarattire on the effects of roughness concentration, 
size variations and relative positions of sediment on ini- 
tiation of motion has been carried out in the following 
paragraphs. 
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Hlffect o f concentration i 

The initiatlonof motion for the bed of maximum rough- 
ness concentration ^ densest bedj has been analysed by several 

investigators such as shields ^.1936^, White ^1940} hurihara 

0 - 

v 1948;, I\;agki \1S^36; and Egiazaraff UsS^J Initiation of 
motion of a single sediment grain resting on smooth bed has 
been analysed by Andreas Hullerj Albert Gyr, & Thonistocles 
Dracos (1971J They studied the forces on single grain in a 
turbulent flow. It may be observed that the effect of 
roughness con. cent ration between the above extreme conditions 
has hot been investigated. Here an attemjit has been made to 
determine the effect of roughness concentration on the state 
of inttietion of sediment motion* 

b . imfect of size vari at^ijm; 

The problem of initiatioh of motion on densest bed of 
roughness eLoments ■ f nonuniform size was studied by Kramer 
■ ^Pef 8 >, UkSWES S^Ref 8>, Cnang i.Ref 8) Indri(Ref.8) ,Aki(Ref.8) 
(Ref. 8) Sakai(Ref,8) ? Satr; and inrihara i^Garde and Raju 
1977.1. The study of ia-'amer, UGtjiS and Cnang indicates that 
the nonuiformity in size Cjiuses an increase in critical shear 
stress* However, the data of Aki and Sato shows that critical 
shearstress for non-uniform sediment is less than the uniform 
sediment. In these investigations the tlireshold condition 
is always refered to the medium size of sediment. The 
initiation of motion of a particular size of sediment in a 
nonuniform bed was studied by Egiazaraff 11967)'. 



It may be noted that the effect of size variation of a single 
grain, subjected to initiation of motion in a given bed of 
uniform size of roughness element, has not been investigated. 
Here an attonpt has been made to study this problem. 

Effect of relative position ! 

The effect of relative position ( degree of exposure) 
of individual roughness element on initiation of motion of 
cohesionless grain was cansidered by Fenton and Abbott. 

<.1977). They conducted the experiments to measure the dimen- 
tionless critical shear stress and its dependence on grain 
protrusioii ' which was found to be very marked. The critical 
shear stress of a grain resting on the top of an otherwise 
flat bed in turbulent stream was measured and found to be 
0.01 which is considerably less than previously reported 
values of 0.03-0.06 for a bed where all the grains were 
at' the same level. 

It may oe noted from the above that initiation of 
motion of a roughness element depends upon the bed charac- 
teristics like relative concentration, relative posit_on of 
roughness elements. These characteristics aodliy the Velo- 
city distr.Lbution near the bed. hence it is of import. •■nee 
to study tne effect of bed characteristics in the velocity 
distribution near the bed. 



The initiation of motion is controlled by the flow 
conditions and the lat«r are aontrolled by the roughness 
characteristics of the bed surface and the state of flow. 

So a brief review of the rouglmess char ectari sties of the 
bed are being carried, out here. 

Wikuradse (1933) found that the velocity distribution 
on a bed having uniform sand grains densly packed can be 
represented in termes of sand grain size K 



in y ^ B ( ) 


( 1 ) 


Where u is mean velocity, Shear velocity and B is rough- 

KV 

ness constant function of ( and velocity destribution 

? V 

near the wall is logarithemic whic±L is universal in character 
and is determind only from wall conditions. 

To characterise the roughness from other type of ro- 
ughnesses ■with different roughness concentration, Schlich- 
ting used equivalent sand grain rouglmess and wall 
Law equation iU can be written as below: 


H = -1 + B ( J .. (2) 

V.J. K —r . — 

* s ^ 

Schlichting, Oloughin and Mac^-o id (1964), David 
(1980) and Sarin (l980) observed that maximum roughness 
occurs not at maximum roughness concentration on the bed 


but at a considerably small value. 
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e- Vartical Shift in the Law of w?i]1 ; 

Hama (l954j expresed the effect of rough wall region 
as a vertical shift from smooth wall log law, and can be 
expressed as bel(.)w : 

^ S " = 1A In 

^ ^ V ® 

or Au 1 In V.*d + ( B - 3 ' ) 

V^= - s r 

^ -1 V^d \ 

Batterraan (1965 J, David, Sarin showed that ( ^ ^ 

K V 

is a function of roughness concentration x ,variation of 
( Y *“ J. ’* } with \ shows a dMinite relationship. 

IC V 

-^^pparent shift in the origin for rough boundrev ( e ) below 
the top of roughness «ilaiient or Theoritlcal ked level . 

The rough wall law behaves as if its origin is loca- 
ted some distance e below t!ie crest of element (More 1951) 
The shift in origin by distance ( ^ ) gives a logarithmic 
distribution near the wall. Perry, Schefield and Joubert 
(1969) suggested a trial and error method for finding £ 

They obtained a relation for closely placed two dimention- 
al strip roughness = 1 ^ 6 

K V ' 

The constant in the equation depends upon the geixietry of 
roughness elanent and has an average v^iLue of 0.40. David 
and Sarin snowed that the constait is a fianction of roughness 
concentration, 
g. Present Investigation i 

The experimental programme is planned to study the 
effect of relative concentration, size and poattioa on the 



ini'Cia'fcion of a single pai‘'bic..e on a bed of unifonu size glass 
beads as roughness elements. Inis study will be based on flow 
paramexers deduced from the mean velocity profiles measured fpr 
the critical condition of a particle. 
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CHAPTER - 2 

EXPERIMENTAL PROGRAMME MD DATA REDUGTION 

The details of the present inTestigation, such as 
description of experimental setup, details of glass bead 
beds, measurement of mean velocity and detection critical 
condition for the glass bead motion are given under the 
experimental programme. 

Deductions of wall law parameters and computation 
of parameters for the initiation of motion have been presented 
under data reduction, 

2,1 Experimental Programme : 

a. Details of Experimental Setup; 

In the present investigation the rectangular flume 
of 15 cm width and 3.0 m length has been used. The system of 
flow straightners are located at upstream end of the flme. 

The water is pumped by a motor through vent urimente red pipe 
which enters into the test section via flow straightners 
for maintaining more or less constant depth of flow vertical 
steel rods of 10 mm diameter and 25 mm, centre to centre are 
provided at down stream end of flume. 

To control the velocity and depth of flow keeping 
discharge constant a drilled Alluminium plate of 15x30 cm 
section used at the discharging end. The discharge is contro- 
lled by sluice valve fitted in the supply pipe. The flume body 
is resting on the lever at centre and is simply supported at 
the upstream end. The downstream end is cantilevered, 
which has been made simply supported, to remove the defect 



cf bending by self load of flume and the live load of water, 
jjoint gauge has been u^ed for measuring the depth of flow 
which can move along the flume length. The test section has 
been chosen at 2m. d/s from the inlet. . A schematic layout 
is shovm in Fig. 1, and plate-2. 

Details of glass beads bed s : 

A 3.Qm long and 15cm wide Aluminium plate has been 
used for sticking the glass beads of 6.5 imn diameter as 
a roughness element to provide the rough bed. A gride of 
lanxiaa. has been drawn for sticking glass beads in a 
regular pattern and getting bed of different rouglmess 
concentration. Five beds of different concentrations have 
been made by sticking glass beads in different spacings 
and arrangraents as given in table 1. For the first bed 
4 cm centre to centre staggered pattern, for the second 
bed 2cm center to center staggred pattern, for the thir^ 
ued 1a.i center to center symmetrical pattern^ for the 
fourth bed 1cm. center to center staggered pattern and 
for the fifth bed densely pa died arrangonent of glass 
beads have been used. After stidcing the glass beads on 
plate, this bed is placed on the bed of flume and is 
tightened oy nut and bolt. Prepared glass bead beds are 
shwwn in Plate* 3 
H easurment of mean velocity : 

A Pitot tube with gauge arrangement was used for 
measuring the velocity profile. A 30*^ inclioied ‘ mcncmeter 
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was used, for measuring the velocity heads. For locatior^ 
of zero position of the probe it is gradually lowered down 
to the position where the pitot tube (probe) is just in touch 
with glass bead top, This is taken as the zero position of the 
probe. The probe has been moved u^werd gradually from this 
position with closer interval for measuring velocity distri- 
tion. 

Discharge measurment : 

A calibrated venturimeter wnich is fitted in supply 
pipe, was used for measuring discharge. 

Slope measurment : 

Trie simple support at u/S end of flume also works as 
slope controlling device as shoim in the Fig 1 
For measuring the slope, the d/s end of flume is shuttered 
by a plate and some water is pumped in to the section, 

WhiHe no flow is there, depth of water was measured at each 
25cm, interval along the length of the channel. From these 
measuments slope of bed was calculated. 
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f. Detection of critical condition for the f^lass bead motion : 

Usually Initiation of motion is detected by the plot 
of Tjf against sediment disch*;rge \,Paintal, 1971). However 
in the present investigation, the approach based on the motion 
of individual rougiiness element has been used. In this approach 
the critical condition is said to have reached when the fluid 
dynamic force acting on the grain increased slightly, the 
grain will be subjected to the motion. The mean velocity 
profiles justat the up stream of rouglmess eltsment are 
measured for this flow condition. The bed shear stress, 
rfmghness characteristics, Ks, and c are computed from these 
velocity profiles. ■ 

2.2 jATA REDUC TIOIT ; 

Velocity measurement were carried out on 4 <iifferent 
densities of glass bead beds, at center line of flume for 
different situations of initiation bf gls-ss b'ca.d mation, Tnese 
velocity measurements are used in deduction of theoretical 
bed level e ) , shear velocity , and roughness func- 

tion like ^is. 

d 

Using the computed values of V^ sheildscritical shear 

PV| 

stress pararaet.^r = and grain shear Reynolds 

■ ^ (p -p) gd 

number R.^^. == was calculated. 

V 

a. Deduction of wall Law par ameters: 

Boundry layerflow can be broadly devided in to two 
regions, namely wall law region and defect law region. 


J7 


The wall law region extended from wall to a distance of 

Y ~ , 

— ~ 0.2 and the remaining portion of a boundry layer is 

repi'esented as defect law region. A portion- of the wall 
law region is represented by a lagarithmic velocity distribu- 
tion which ranges from y/5s:0.05 to 0.20. Lagarithmic velocity 
cistribution is universal in character snd it is the function 
of wall characteristics like bed shear velocity roughness 
function and theoritical bed level e , These characte- 

ristics are detemined from logarL thixic ve3.ocity distribu- 
tion near the wall. 

b, . Theo ritiG?J. . Bed l evel ( ^ ) 

In the case .^f rough wall flow, it has, been found 
by iioore Perry, Joubert k 1963 ) and Liu etal 1^1966} 

that the logarithmic profile still exists provided that the 

Y origion is located at; a distance ^ below the crest of the 

roughness element i.e. Y^s= where Y is measured from the 

top of rougimess eleraent. In the present analysis^ a graphical 
method is used for the determinatitm of e and shear velocity 

The basic steps required to determine £ and using 
the graphical method are; 

1, The orgin for Y^= Y+'c is chosen by trial and 
error to give the closest approximation to linearity in 
the seiai lc;garithmic graph. 

2. The best straight line is drawn through the 


resulting ' plot| 
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^ahear Ve loci ty V,^ . 

Shear velocity is defined as =/ .Jl. in i/y^'hich 

P 

IS bed shear stress and p is fluid density. The slope of 
the best fitted straight line in the plot of u against In 
^t the value of The value of V,_^ is computed for 

the Larman constant k = 0.40 from, above computed slope. 

The bed shear stress is computed using the relation. 

/_Lin vhich p is density of fluid. 

RougjQnjess funct ion 

The rougliness effect represents the doivn word 

displacement of the logarithmic portion of the plot of rough- 
wall law u with \Y' ¥. c ) from the smooth wali log law. 


The smooth wall log law is represented by equation as 


below. 


u 1 YV. ^ . 

7 ~ ^ 5.5 


And wall low for rough bed may oe represented by 


relation. 


i In LilL.L-U.-* .^5,5 


The vertical shift of rougii wall l(jw eqh (2i from 

smootn wall law eqa (.1) is measured from the plot u with 

% 

c ) j_n Fig 4.A . 

V 

Eq uiva len t sand grai n Roughne ss <,Ksj 

For hydrodynamically smooth and rough boundaries, well 
known velocity distribution equation given by Prandtl and 
Karaman as wall lawi 

V 

For Smooth boundries lU , = 5.75 log 5.5. .. .aI,) 

V* V 


^ ^ V 

For rough boundries: =5 75 ^ ®*5 


...( 3 ) 










2i 


From these two equations UJ & (3) we cen i- et 

3^.0 


Au 

V'- 


u , 

'Vis 


V'^ks 

5.75 log V 

ks can be calculated from this equation using the 


./M) 

kid 




values of au and V 


Parameters for the in itiation of motion 

initiation of motion of sediment was represented by 
Shields U936^ in terras of shear stress parameters 

f>l«' «<r 

* w - — and grain shear f^eynolds number 

k— -1^ Sd 

R.«. =V.ft d . using the computed value of V^, and measured 

V 

values of grain diameter d, and density of ' glass bead 
and known value of p for water at' given 
rature, the values of and t have been computed. 


Ejcperimental data is given in Appendix, 
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CHAPTER-^ 


AI^TAI^Y SIS op REg[[L'rs 

The experirflontal results have been analysed in the following 
sequence; 

1, Effect of rougliness concentration on the mean 
velocity distribution and their roughness scales, 

2, Jiffect of Roughness concentration on the initiation 
of motion of glass bead; 

3, Effect of etze variation of single part i cl e( glass 
bead) in given bed on its initiation of motion, 

4, jiffect of relative position of single particle on 
its initiation of motion, 

1.1 Mean velocity D ist ributi on; 

The mean velocity profiles measured near the bed 
for the different roughness concentrations are presented 
in two methods as indicated in the introduction. In the 


first method, the velocity profile is represented interms of 


shift in velocity namely Au , as 

f In 5.5 - 'S 

■# 

Incorporating the velocity shift Au with velocity 


profiles as 

u+Aii In ^ 5.5 ...Ua> 

and plotting yi Au against kIZL). ^ in fig. 5 for all beds 

V , V 

of various roughness concentrations, it may be observed 


that the experimental data agrees fairly i^ith Eqn, la. 

In the second method of presentation, the velocity 
profile is represented interms Nikuradse's equivalent sand 
grain rouglmess as 







i 2 ) 


2U 


u I . 

h‘ 


Ihe experiment 

T 
k 


'iata plotted in the form of u against 

i-f. u • ^ . 

— = V’^ e •^^S* t> agrees fairly with the eqn (2) 
k •, ~ 


'■s 
1.2 


RoUKiines s concentration on 
lioughness sca l^ 


The roughness scales used in representing the 

velocity profiles as shown in fig 5,6 are — ~ and . 

V* d 

These scales are found to h® th® function of zoughness concen- 
tration. The roughness scale plotted in the form In civ..^) 

it V 

against roughness concentration x as sliown in Fig. 7 is 
found to increase with intarease of x upto x =0.26, and then 
decreases with further increase in x . The data of gciilichting 
Oy37i David 0980) and Sarin ^1980) are also incorpated in 
7» which showgi fair agreement with the present investi- 
gation. David present an explsoation based on the variation 
of effective area of rouglmess eleraent for this particular 
characteristxc variation. The roughness scale presented interras 

of k/d is plotted against X in the Fig. 8 alongwith the data 

s 

of Schlichting (1937)» David v1980) and Sarin (l980). It is 

K 

seen that tiie value increases with increase in X , attains 

d 

a maxiraim value around X = 0.26 and decreases further with 

increse in x . This characteristic relrjti n of k /d with 

s 

X is similar to the variation of ( ^ - 7 In ^~*) with X 
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Theuriticpg B.^d 1 e 

The app=irent origin or the thaoritical bed level 6,f 
the logari-thLUc val.ocity distribution situated below the top 
uf the roughness element xs of impartance in the computation 
of fluid dynamic force acting on roughness eLenient, The 
location of this origin depends upcm the size and roughness 
concentration of gedement bed- B^'inco and PartheniaKies 
0971^ gave, a relation e= 0.25 ^sffor the location of apparent 
origin below the tup of bed surface in which sediment grains 
are densely packed. The Theoretical bed level found to vary 
with roughness concentration on the' bed, as sho\\?n in Fig. 9 
It may be observed that as the rouglmess co>ncentratiun decreases 
the thearetical bed level moves "feo words bottom of roughness 
element-. Data of David and Sarin have also been in corporated 
in this' figure. In the same figure, the surface level of the 
rough bedhe§ «lso been plotted and it may be noted, thft 
the theoritical led level a>incides with the surface level 
of the bed. 

c. General observations 

The roughness concentrati ,>n has a major • influence 

K 

on the roughness factors. The r<;ughness factor is the siua of 

d 

the rouglmess contribution by the rouglmess elements and smooth 
portion of the bed. The fv,n.ctioncHl form of roughness contribu - 
tion by the floor and roughness elcanentg can be written according 
to Einstien (1950J as fellows f^CTutal ^ea of bedj = 

Ahrea of smooth portion of bed}'+‘ lArea of Roughness element) 
in which T is C'jmbined shear stress of smooth portion of bed 

^ O ' ■ ■■ 

and rouglmess elements on the bed. Suffix 'f and 'r' represent 
smooth and rough bed respectively. Dividing the above expression 
by totalarea of bed , one can re-write it in terms of concentration 
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0 \)f ^ ^ (l) 

Since the projected area of roughness elment is 
equal to Klan area of roughness element and it can be related 
as Xf = 1 - Xp .usi.ng this relation in equation ( 2 ), and 


denotiiqg = x 


^Or To 


= J_f 1 

To.. ^ X ~ X 


(5) 


— remains fairly constant fur given flow condition . 

° To 

The ratio ^ will he a function of x and this functional 
0 

form has been plotted in fig. 10. It maybe ob^rved that t^^ 
increases with decrease of represents the shear stress 

ofiered by single roughness elanent.. may be considered as its 
drag f 01 ce. This may be related with the drag coeffecient 
and velocity of flow Ug.at the top of rougJiness element » as, 
given belowi 


Drag^ c'Force _ 1 pu ^ 


...U; 


Or projected area. 2 "'C 

The velocity u^ may oe computed by using wall law equation. 


u 






+ B) 


...(5; 


knowing the values of e , ks and = 8.3 for rough turbulent 
fLo'vs, the Ug , velocity on the top of roughness element is 
computed. C^^Drag coeffecient was computed by using this 
velocity UgSnd has oeen plotted against X in Fig. 11. %e 
data of David ^1980.), Sarin (l980) for glass beads and David 
for sands has also been in corpora-ted in this figure. The 
data of Schliciiting ^1937) has also been incorporated. Now it 
is observed that the data of the present investigation agre 
With David and Sarin expirimental data. Further it is' noted 
that C£) remains fairly constant for i < 0. 10 and decreases 
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with furtiier increase xn x. .ji ernperical equation fer 
this curve may be written as vl- Xwnere G^- 

0.45-) . The v, lue of for X <0.1 fairly coincides with 
the values calculated, by Weil L. Coleman ^,196?^ in the 
rough turbulent zone. The computed value of from Schli- 
chting data do not agree with the present investigation. 
Reason may be attributed to the absgence of information on 
theoritical bed level and hence on the velocity Ug . 

The variation of drag co-efiecient with roughness 
conceiitration X infer the variation in the magnitude of 
fluid dynamic force acting on tlie rouglmess element. The 
fluid dynamic force may be resolved into drag force parallel 
to flow and lift force perpendicular to flow-Chepil points 
’hat the ratio of the magnitude of lift force and drag force 
for rouglmess element resting on a bed remains fairly con- 
stant. Hence it may be stated that tlie fluid dynamic 
force acting on the roughness elieiient decreases as its 
concentration on the bed increases. This may result in 
the need of more fluid velocity near the bed, or the 
shear stress to cause the grain t-mmove as the grain a)n- 
centration increases. 

2 . Eff ecct _of _RougImess conce ntration on Initiation 

o f motion : 

Shields proposed that the dimensionless shear stress 
T.^ = — 1 ) necessary to move the grain on 

stream bed dependes on the grain shear Reynolds number 

However, it has been found by White v1940) that 

•■fr- V 



fact..;r of relative packings of roughness elements alter the 
initiation criterion proposed by Shields (1936; 

The initiation of motion of 6.5mm diameter glass 
bead for difierent density beds ^= 0.0414, 0.166, 0.332 
and 0.753 is ooserved. The experimente were repeated for 
number of times by varying the slope, discharge and depth of 
flow to obtain the critical condition. The velocity distri- 
Dution were carried out for the condition when the roughness 
element moved out from its position. From this v^ocity 
distribution shear velocity V ^ and other roughness parameters 
were computed. Using these computed values the Sliields 
parameter t,. and are computed and plotted in Fig 12. 

Ihe data of David (l9B0j and Sarin for glass bead and sand have 
also been incorporated* It can be observed from this figjire 
that all data of varying density, falls systamiatically in 
between the data of daisest and smooth bed, and it is observed, e 
that as the rougimess density increases the curve for isi.iootii bee 
movesto^^ords the curve for densest bed. The data of densest bee 
^ = 0.753 of present investigation falls very close to 
Shields curve. The shear stress offered by single grain 
has been ploted against in Fig 13 and it can be 
observed that data for all roughness beds falls in the 
limit of Shields curve. 

5. Ef fect of si :e variation of single particle (gl ass 

bead) in given bed o n its initiation of motion . ■ 

To study the efiect of size variation of single roughness 
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element in a given rough bed, on initiation of nntion, the size of 
sediment (d^) subjected to initiation is changed and for this state 
of initiation of motion, shear stress is calculated and plotted 
against in Pig 14 along with Shields curve. The data of David 
for different rough-beds have also been incorporated. It can be 
observed that as increases decreases and the curve shifts 

towards Shields curve as the roughness density of bed (^) increases. 
The shear stress offered by single grain t is calculated and plotted 
against in Jig 15. It can be observed that 

T 0 = -1 

d^ 

Calculating the parameter ^ v.,,- 

'^TX u X 

in lig l6 , It can be observed from this plot that the value of this 
parameter is almost constant and very close to Shield q^rve in turbulent 
zone for present data. 

4 . Effect of relative protrusion of single particle on its initiation 
of motion 

In order to study this effect, roughness element vtoich is 

subjected to study the initiation of motion is placed at different 

levels P = d/4, d/3, d/2 above the rough surface, and condition of 
s 

initiation of motion is determined for these positions on the bed of 
roughness concentration X = 0.753. 

The shear states \ is calculated for this condition and plotted 
against j^/d in ilg 17 and the data Poaton and ibWtt has also been 
incorporated, it oan be observed that present expeiimaatal data agrees 
fairly with Penton and Abbott curve. 





shields curve 
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For different Pg/d ratios, t* with R# is plotted in Pig. 18 
it may be observed that as Pg/d increases t* values 
decreases, 

Disoussiong * 

Prom the analysis of experimental data, it 
may be stated that the value of at the threshold 
condition decreases from the Shield’s critical value, 
with decrease in particle concentration (roughness 
concentration) , with increase in particle size, and increase 
in in relative position of a particle with respect to 
general bed level. An attempt has been made to explain 
this characteristic variation using equilibrium of forces 
acting on a particle. 

Prom the equilibrium of forces acting on a glass 
bead, for its initiation of motion, it may be shown that 

Drag force > (submerged weight-lift force) tan 0 
where 0 is angle of response 

2 2 

Drag force = 1/2p Ug d 

Lift force = 1/2P Ug 0^ d^ 

3 

ubmerged weight =(Yg'-^) G 2 d 

Substituting these expressions in equilibrium equation 
after rearranging, one gets 

l/2p u| (Cx>-Gl tan 0 ) > (Yg~Y) 02 ^^ tan (6) 
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Ihe velocity Ug at tlae top of the grain may be 
computed from equation 

Ip 

Ug = V, ( _ In ^ + B ) 

Value of constants and depends upon the 
shape of the particle and for spherical glass bead 
0 ^ = te/4 and O 2 ~ 7 c/ 6 . The angle of repose 0 depends 
upon size and angularity of the particle. For spharical 
glass bead of 6,5 mm diameter 0 is taken as 45 °. 

For initiation of motion of glass bead, the 

threshold condition may be written as 

2 

V* 1.335 

- .... ; > — p 

(Yg -Y)d ( 5.75 log ^ +B)2(0jj-0p) 

s 

The threshold condition depends on the position of rough- 
ness element with respect to theoretical bed level ? 0 > 
the roughness characteristics of bed Kg, the roughness 
constant B, and coefficient of Drag Op and coefficient of 
Lift Cp of the particle. 

These parameters vary with the roughness concen- 
tration, size of the particle and relative position and 
the state of flow. 

As the roughness concentration decreases from the 
densest case, the value of Drag coefficient Gp and the 
relative position of theoretical bed level e/d increases 



4-1 (b) 


shown in Pig. 11 and Pig. J respectively. The value of 
roughness scale IC.g/d depends upon the roughness concen- 
tration as shown in Pig. 8. Variation of lift 

coefficient is not known. However, making assumption 
that the variation of lift coefficient is similar to 
Drag coefficient, it may be shown from Equation 7 that 
T* decreases with decrease in roughness concentration. 

This relationship has been observed experimentally in Pig, 12, 
. When the size of particle (di) subjected to 

threshold condition, incomparision to the size of roughness 
eloments(d} increases, the value of P0/d also increases. 

This effect causes the decrease in t* with increase in 
dj^/d, This variation has also been obsert^d experimentally 
as shown in the Pig. 14 . 

The increase in relative position (elevation) of 
the particles, causes increase in P@/d, This results in 
decrease of t* as P0/d increases in equation 7 which can 
be observed also from Pig, 17 . 

The effect of decrease in roughness concentration, 
increase in size and increase in relative position causes 
the threshold values to decrease from the threshold value 
given by Shields for the densest bed. By representing 
the threshold condition based on the shear stress offered 
by Particle subjected to motion (t* ) instead of shear 
stress of the bed x*, it can be seen from Pig. 13 that 
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values a^p?ees with. Shields values fairly well for 
all roughness concentrations. This kind of analysis could 
not be shown for the variation of relative size and 
relative position of particles because. the value of 0 ^ 
was not calculated for the above cases. Prom this study 
it • can be brought out that the knowledge on the variation 
of Oj), , e and Kg are needed to predict the threshold 
condition of a particular particle. In this direction, 
attempt should be made for determining 0^ and of a 
particle resting on bed for different relative positions 
and for different relative sizes. 



CHAPTER^ 4 

cbi?ax..U:jioFa 




A study on the effect of roughness concentration, 
size variation and relative position of a single particle on 
the initiation of motion has been made. Roughness concentra- 
tions A = 0.041, 0,166, 0.332, 0.646 and 0.753 are us'^d in 
this study, fho size of particle 'vdj_J in comoarision of rough 
ness element size (d/, namely d^/d was varied. The relative 
position of particle measured from the theoretical bed level, 
namely P/d was varied. 

C' 

The effect of size and relative position of a particle 
on the threshold condition are carried out in the densest 
roughness concentration of the bed. For each threshold condi- 
tion, the mean velocity was measured at the test position. 

All thc' rel event rougliness parametors are de^u«fc, 0 d from 
these mean velocity distributions. Based on these experimen- 
tal observations, the following amclusions are drawn: 

\ b ) C^j-ofi'ecient' of drag; The Co-effecient of drag 
of ri)Ugimoss eltanent increases with decrease in rougliness 
concentration and remains fairly constant for X ^0.1. 

\,b) Theoretical bed level: Theoretical bed level ' e ' 
shifts from the wall surface on which roughness elements 
are attached towards the top surface of roughness elements 
as roughness concentration increases. It was found that the 
theoretical bod level coincide with average surface level 
of a rough bed. 





i,c; Initiation of motion; 

\i; Thrssiwld condition represented by Shields 
■parameter for rough bed condition decreases with dec- 
rease in r!.)Ughness concentration and attains a value of 
= 0.01 for the case of single roughness element on 
siiiooth bed. liy considering the shear stress offered by 
individual roughness elemmtsj it is jiiowi that the critical 
shear stress coinside fairly with Shields curve. 

viij It is obser’eved that the shear stress of a single 
rougi'uiess clement ' t^resting on the bod, decreases with 
Lncreaso in size of roughness elanent. The effect of size 
variation taken in the form remains fairly cons- 

tant and this value agrees ’i\/ith Shields critical shear 
stress value. 

Iviiij The critical shear stress decreases with 

increase in the relative protrusi'jn ]|/d in the rough turbu- 
lent ztjno. Present experimental results agree with the 
Fenton and Abbott results. 
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d 

6.5 mm 

A 

APPEiroiX 

. 041 4 Kg/d 

e/d 

= 0.853 

= 0.770 


Run B 

M 0 . cm 

5 

cm 

V xK/''^ 
cm^/se'j 

Y.-, 

' »\ 

cm/sec. 

^aug 

cm/sec. 

R*x10“^ 

xIO"*"^ 

1 2 

3 

4 

5 

6 

7 

8 

9 10 

1* B,7G 

3.90 

B.80 

2.55 

31.66 

1,88 

0.602 

0.32 

2. 6.76 

3.97 

B.20 

2.26 

29.05 

1.80 

0.548 

0.36 

3. 10.94 

4.43 

B.20 

2.37 

29.36 

1.88 

0.602 

0.28 

4. 12.08 

3.97 

8.20 

2.27 

28.63 

1.88 

0.602 

0,22 

5. 10.40 

4.47 

B.10 

2.23 

29.54 

1.87 

0.584 

0.29 


X = 0.166 Kg/d = 3*176 

e/d = 0.61 


1 2 

3 

4 

5 

6 

7 

8 

9 10 

6. 11.28 

4.92 

B.84 

3.35 

31.17 

2,59 

1.20 

0. 30 

7. 8,88 

5.50 

B.B4 

3.32 


2.57 

1,18 


8. 12.41 

4*42 

8.80 

3.50 

33.70 

2.59 

1.31 

0.31 

9 . 9.34 

5.42 

8.80 

3.39 

33.78 

2.50 

1 . 23 

0,34 

10. 9.35 

6.40 

8,70 

3.83 

hmf 

2.89 

1,57 

- 

11. 11.54 



5.42 

8.70 

3.39 

31.21 

2.72 

1.23 

0,27 
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0.552 KgA = 5.37 

e/d = 0.46 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 2 

9.92 

6.3 

2 1 80 

4.62 

49.18 

3.41 

2.29 

0.50 


15 

11.75 

6.3 

9,30 

4.33 

39.91 

3.06 

2,01 

0.37 


u 

B.45 

7.3 

9.20 

4.39 

41.73 

3.10 

2.07 

0.46 


15 

11.75 

7,3 

9.20 

4.28 

40.89 

3.02 

1,96 

0.38 


16 

10.36 

7.3 

10.50 

4.35 

43.80 

2*68 

2.01 

0.43 


17 

10.95 

7.3 

10.50 

4.50 

44. 24 

2.68 

2.01 

0.43 


18 

7.47 

6.3 

10.20 

4.33 

42.73 

2.76 

2.01 

0.50 


11) 

9.25 

7.3 

10.20 

4.33 

47.15 

2,76 

2.01 

0,50 


A » 

Alwnlnlura 

gr= 

/?ravel 


P3/d 

= 0.25 



S 

,31 ass 


P = 

Plastic 

Kg/d 

= 1.7 







X =:0. 

, 66 

e/d 

= 0.31 



i... 

2 

3 

4 

5 

6 

7 

8 

9 

10 

20 

6.55 

6,55 

3.10 

5.94 

78.93 

4.77 

3.79 


0.65 g 

21 

7.50 

7.50 7.90 

3.65 

66.43 

2.94 

4.04 

0.77 

0.625P 

22 

7.86 

6,10 

7,80 

5.83 

63.68 

4.68 

3.65 

0.73 

0.650g 

25 

6.5B 

6,58 

8.00 

6.27 

78.73 

2.90 

3.73 

0.9s 

0.750gr 

24 

5.45 

5.10 7.70 

6.60 

84.84 

6.00 

3.73 

1.16 

O. 7 OOAI 
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X = 0.66 




1 

P. 

3 

4 

5 

6 


8* 


io" 

25 

0.47 

6.15 

S.O 

4.11 

48.57 

5.58 

2.09 

0.50 

0.65 g 

26 

9.17 

6,15 

7.9 

4.37 

49.30 

4.17 

1.82 

0,53 

0.75 gr . 

27 

8. 92 

B .92 

8.0 

2.77 

43.75 

2,20 

2.34 

0,47 

0.625 P 

2 B 

11. 1 H 

6.90 

H.O 

4.37 

32.80 

3.55 

2.05 

0.31 

0.650 g 

29 

1P.90 

1 1,00 

7.7 

2.73 

33.71 

2.26 

2.27 

0.30 

0.625 P 

30 

10.7P 

. 00 

7.9 

4.59 

38.31 

4.43 

2.01 

0.37 

0.750 gr 

31 

H. 25 

7.95 

7 , 8 

4.92 

60.55 

4,31 

2,07 

0.67 

0.70 A 1 






X s 

0.66 









Pg/d 

=0.75 




1 

2 

5 

4 

5 

6 

7 

8 

9 

10 

32 

10.50 

7.00 

B.OO 

3.28 

49.61 

5.09 

1.02 

^ T 

0.48 

0.75 g 

33 

9,55 

B ,20 

8,00 

3,43 

50.51 

3.00 

1.01 

0,52 

0.70 A 1 

34 

9,73 

9.10 

■% o 

0.00 

3.06 

46.62 

■2.4§ 

1.01 

0.48 

0.65 g 

35 

11,50 

7.05 

8.00 

2.11 

(Td* 

28 l 92 

1.68 

1.36 

0.27 

0.625 P 
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l = 0.753 

Pg / d - 0.00 


1 

o 

3 

4 

5 

6 

7 

8 

9 , 

10 

H 

6,11 

5.35 

8.00 

6,42 

63.12 

5.22 

4,41 

0.82 

0.65 g 

37 

6 , OK ^ 

5.85 

8.00 

6.59 

52.81 

5.37 

4.12 

0.68 

0.75 gr 

38 

5.44 

4.85 

8.00 

6,93 

63.09 

6*06 

4.12 

0,86 

0.70 A1 





X 

» 0.753 








Pg / d « 0.65 

cm 




1 

2 

3 

4 

5 

6 

7 

8 

9 , 

1 0 

39 

9.48 

9.48 

8.00 

4*48 

2,30 

3.46 

0,405 

0,62 


40 

9.58 

6.90 

8.00 

5.03 

3,02 

3,67 

0,48 

0,74 


41 

10,95 

7.10 

7,90 

5.76 

3.91 

4.36 

0,53 

o 

CD 


42 

6.11 

5.35 

8.00 

6,42 

5,21 

4,41 

0,65 

1,00 


43 

4.30 

4.30 

8,00 

6.93 

6,32 

4,18 

0,80 

1,23 


4 4 

. 4.22 

4.22 

7.95 

10.92 

13.48 

8,31 

1,00 

1.54 


45 

4.15 

4.15 

7.80 

12.06 

15.46 

10.14 

1.10 

1.69 
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kymmw . vmvues of RoufflNEss functions md shear stress 

for (1LA.3S BEAU BED 


Unifom siTie d == 6,5 nm Arranged in a Systamatic pattern 


SI. 

No. 


X 


e/d 




AU T*d 



1 0.0414 f ’^.77 

2 0.1660 0.61 

3 0 .' 5-520 0.46 

4 0,6640 0,51 

5 0.7530 0.23 


9.7 

0.853 

189.56 

13.75 

3.176 

257.71 

14.48 

3.750 

292.42 

14.10 

1.700 

486,01 

13.30 

1.268 

539.32 


6.974x10“^ -3.39lf 

1.2462x10"^ -0.114 

2.150x10*"^ +0.300 

3.648x10"'^ -2.500 

4.220x10"^ 


-2.400 


